Stochastic theory of quantum vortex on a sphere 
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A stochastic theory is presented for a quantum vortex that is expected to occur in superfiuids 
coated on two dimensional sphere S^. The starting point is the canonical equation of motion ( 
the Kirchhoff equation) for a point vortex, which is derived using the time-dependent Landau- 
Ginzburg theory. The vortex equation, which is equivalent to the spin equation, turns out to be 
the Langevin equation, from which the Fokker-Planck equation is obtained by using the functional 
integral technique. The Fokker-Planck equation is solved for several typical cases of the vortex 
motion by noting the specific form of pinning potential. An extension to the non-spherical vortices 
is briefiy discussed for the case of the vortex on plane and pseudo-sphere. 

PACS numbers: 05.40.-a, 67.25.dk, 67.25.dr 
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I. INTRODUCTION 

The study of vortex dynamics has been one of central 
subjects in classical physics[l|. The modern interest of 
vortex phenomena is motivated by the quantum fluids, 
namely, the quantum vortex in various contexts, which 
are considered to be defects inherent in the complex order 
parameters describing the quantum fluids. Specifically, 
the dynamics of point vortex has been extensively inves- 
tigated so far [2-0L, The dynamics of vortex string has 
been also studied Q . Besides the conventional vortex in 
superfluid He4, the more complicated textural structure 
in He3 has been explored ( @ ) . Furthermore even the vor- 
tex in the cosmological scale has been explored(e.g.[7]) . 

In conventional treatment the vortex is an object de- 
fined on Euclidean space of 2 and/or 3 dimension. The 
vortex defined on curved space rather than fiat Euclidean 
space is possible to occur. Indeed, the classical counter- 
part of non-planer vortex has been given in several con- 
texts starting with the purely classical hydrodynamical 
procedure, see e.g. [l, (oj. 

On the other hand, it is known that in condensed mat- 
ter physics the vortex motion is influenced by random- 
ness caused by e.g. the temperature fluctuation and the 
presence of the interaction with the various sorts of im- 
purities. For such cases, we need to treat it within a 
framework of stochastic theory. Indeed, the stochastic 
theory of vortex has been investigated in connection with 
the superconductivity • 

The purpose of this note is to present a general for- 
malism for stochastic theory of quantum vortex that is 
constrained on the two dimensional manifold; specifl- 
cally, two dimensional sphere S^. The essential point 
is to follow an analogy between vortex and ferromag- 
netic s pin : that is realized for a single magnetic "do- 
main" In order to achieve this, we flrst derive the 
equation of motion for the vortex center on these man- 
ifold; the canonical (Kirchhoff) equation by using the 
time-dependent Landau-Ginzburg theory 12], that is de- 
scribed by the complex order parameter on (Section 
2) . Then the equation of motion for the vortex is con- 



verted to the Langevin equation by noting the analogy 
between the spherical vortex and the spin in the presence 
of the random force (Section 3). From this analogy, the 
Fokker-Planck equation is derived by adopting the func- 
tional integral technique. The functional integral used 
here is borrowed from the procedure that was developed 
for the statistical theory of wave propagation in random 
media [Tsl. [l^ . In section 4, we show the solutions of 
the FP equation for a typical form of pinning potential. 
In the section 5, we give a brief sketch of an extension 
of the spherical vortex to its variants; the planer vortex 
and then pseudo-sphere vortex. 



II. EQUATION OF MOTION FOR QUANTUM 
VORTEX 



We start with the time-dependent Landau-Ginzburg 
Lagrangian, that is expressed in terms of the order pa- 
rameter tp: 



1^2^^ dt 



c.c.)-n{ip.ij*)]da. (1) 



This form is the same as that was used for the superfluid 
which is deflncd on a flat plane [l^ . Here tp{r,t) is de- 
fined on 2-sphere S"^ with radius |r| = a and the integral 
/ da is taken over S"^, the coordinate of which is denoted 
by r written in terms of the Cartesian form. H represents 
the Hamiltonian density which consists of the kinetic and 
potential energy; 



n = — VV'*VV'- 
2m 



(2) 



where V means the derivative with respect to the spheri- 
cal coordinate and m means the mass of constituent par- 
ticles (e.g. He4 atoms). 

Let us consider the case that a single vortex is created, 
for which the order parameter is written by a form such 
that 



V'(r,t) =VXr-R(i)). 



(3) 
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Here R denotes the time-dependent vortex center satis- 
fying |R| — a. With this parametric form we have 



-^(V/V.V-cc.)^, (4) 



where use is made of the relation: 

together with the relation = — Vr (that is, the the 
nabla with respect to the vortex center coordinate tuns 
out to be the derivative with respect to field argument 
r). Now we adopt the polar form of the order parameter 



ip = y/pcxp[i—a{r - R(i))]. 



(5) 



Here p = p(r — R(t)) is the density, which vanishes at the 
vortex center r = R and tends to constant value: p = po 
outside the coherent length. The first term of L (denoted 
by Lc that is called the canonical term) is calculated as 

j • —da = / {j, -X+jyY+j,- Z)d<j, (6) 

where mass current j = mp\/a = mpv is defined. On 
the other hand, the Hamiltonian term becomes 



H = l[lmpv^ + j-{Vpy+Vip)]dcT 



(7) 



The first term means the fluid kinetic energy and the 
second and third terms are the internal energy, which is 
quoted as U{p). The exphcit form will be given later. 

Now the equation of motion for the vortex center is 
derived from the Euler-Lagrange equation: 

d ,dLs dL 

The contribution from the Hamiltonian is 
dL _ d f _ OH 

dR-mJ 9R- 

On the other hand, the canonical term is manipulated as 
follows: First by differentiating the under integral sym- 
bol, we have 



d ,dLcs 
di^'dX- 



dt 



■da 



together with ^ = + + and 



dY 



dz • 



die f ,dj^ ^ , djy^^ , (9^2 



Hence we obtain 



(8) 



where use is made of ^ = — By evaluating the 
RHS of ([5]) (see Appendix A), we arrive at the Kirchhoff 
equation of motion for the vortex center 



/iri(R X 



dR 

It' 



dH 



(9) 



with 51 = mpQa? and replacing as i7 — > H / a. p means the 
"vortex charge" , which is introduced a multiplied factor 
in velocity field (Appendix A). Eq.® coincides with the 
equation of motion for classical vortex on sphere [1, Q , 
which is derived in the framework of classical fluid me- 
chanics. This equation has a simple dynamical meaning; 
the balance of two types of forces; the left hand side rep- 
resents the Magnus force and the right hand side is the 
gradient force coming from the Hamiltonian (that plays a 
role of potential energy in analogy with classical particle 
mechanics). This equation of motion can be written in 
an alternative form: 



pVl- 



dR 

It 



R X 



dH 



(10) 



The equation of motion is rewritten in a form of Hamil- 
tonian equation of motion. Namely, by using the compo- 
nent for spherical vector basis 

dR • 

— = (0,e,sine$), 

the equation of motion is written in terms of the angular 
form: 

^ne ^ -i^ll, pm = --^H. (ii) 



sin 8 9<i> ' 



sin e de ' 



This is of the same form as the equation of motion for spin 
(see, e.g. [H, HB]). The effective action corresponding to 
the equation of motion may be given as 



5'= [pQil - cos e)^ - H]dt. 



(12) 



Quantization of vortex charge: We examine a special 
feature of the above form of action function; namely, the 
line integral. As a particular case, we consider the in- 
tegral along a closed loop, say, C. This can be written 
by the surface surrounding C; where there is ambiguity 
choosing it, 

/ ^17(1- cos e)d$ ^ pn I sinOded* 
Jc Js 



-pn / sin eded$ (13) 

Js 



which consists of the upper or lower surfaces, respectively, 
say S and S that are complement each other; S+S = S^. 
This ambiguity is related to the gauge choice [l^- Upon 
quantization, the ambiguity is solved; that is expressed 
as an exponentiation; 

exp[^pQ / sin9d9d$] = exp[--^//ri / sin8d6d<l>] 
I*- Js J s 
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then we obtain 



A. Spin analogy and the functional integral 



/ifi / sin QdQ A — 2mTti (n = integer) 



S2 



which leads to the Dirac monopole quantization: 



ATTflil, EE ^^2' 



(14) 



where M = Airfl, which means the mass of fluid on the 
sphere. Equation ([T4l) is just the same as the relation be- 
tween electric and magnetic charge for the case of quan- 
tized monopole. 

Hamiltonian term; pinning potential: 

The Hamiltonian term is evaluated as follows: First we 
note that the fluid kinetic energy together with the 
second term which is given by the gradient of the density 
profile. These two are independent of the position of 
vortex center for the case of single vortex. We need to 
settle a coordinate that is relative to the vortex center. 
Then the term satisfying this criterion is the last term 
V{p). In order to express this explicitly, we consider the 
interaction energy with the pinning centers that is built 
in the superfluid, which is written in a form 



U = / V*(r)V'(r)V"(r)dr 



(15) 



and noting that the density profile is given as |V'(r)P = 
po + p(|r — R|) where po is the uniform background term 
and p describes the vortex profile, which may be given 
by a Gaussian form with the central peak at the vortex 
center. The interaction ^(r) can be chosen such that it 
couples in a contact form; namely, that may be given by 
the delta function: 



V{r) = Vo6ir - a) 



(16) 



with the pinning center a, then we have U — Vqp(II — 
a) up to an additional constant term coming from the 
uniform background. If we take a = (0, 0, a) as the north- 
pole: it turns out to be 



U = l/op(aVl(l-cose|) 



(17) 



where O denotes the angle between the vortex center and 
the north-pole. In the case that there are several pinning 
centers located at afc, U is given by 



N 



U = Y.^ok{pm-ak\) 



III. 



k=l 



THE LANGEVIN AND FOKKER-PLANCK 
EQUATIONS 



In this section we consider the stochastic equation for 
the vortex motion by following the analogy between spin 
equation of motion. 



As is remarked in the previous section, the equation 
of motion ^TU\i is of same form as the spin equation of 
motion; the vector R, plays a role for the spin vector, 
which we call the pseudo-spin. Here we introduce J — 
JR, and J = pQ corresponds to the spin magnitude. In 
terms of J, the above equation is written in a scaled form: 



dJ 

dt 



— = -J X 



dH 



(18) 



The right hand side of the equation represents the torque, 
namely, ^ is nothing but the magnetic field acting on 
the pseudo-spin, which may be called "pseudo-magnetic 
field" . From physical point of view, the analogy between 
spin and vortex is not surprising, since the vortex may 
be regarded as a region where the angular momentum of 
the fluid is concentrated. 

Furthermore, following the well known fact in the spin 
theory, it can be extended to include the effect of dissipa- 
tion: such that VH is replaced by VH+rj^ in P^ : with 
77, the dissipative coefficient. Having solved the equation 
for fl, it follow 
equation" pl|): 



for it follows (which is known as " Landau-Lifschitz 



A(J) = , 



where we use the notational convention: 

A( J) = TT—r J X — - - 77 — - , 



(19) 



(20) 



We now treat the Brownian motion of point vortex 
that is caused by random effects of several origins; the 
temperature fluctuations, and inevitably existing impuri- 
ties. Here we follow an analogy that is expected to occur 
between a single vortex and a ferromagnetic spin which 
can be observed in a single ferromagnetic domain fl]| . 
For the case of ferromagnetic spin, the simplification is 
adopted; the random thermal fluctuations have a corre- 
lation time much shorter than the response time of the 
system. The response time of a single domain particle is 
of the same order as the reciprocal of the gyro-magnetic 
resonance frequency. However, the case of the vortex on 
a sphere, it is not feasible to estimate the ratio between 
two characteristic times, so we adopt it as a working hy- 
pothesis. 

In what follows we consider that the fluctuation effect 
gives rise to the randomness of the pseudo-magnetic field; 
which may be written as b, so we replace ^ — ^ ^ + b. 
Here b(t) is assumed to be the Gaussian white noise as 
a random uncorrelated magnetic field. Hence we have 



d3 



A(J) -c(t) , 



where c{t) the right-hand side, that is 
1 



cit) 



1 + J2 



(J X b - riJ-'h) 



(21) 



(22) 
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which is a combination of random force and torque. As 
a result of of the random uncorrelated features of the 
function b, we assume that c can also be expected to be 
Gaussian white noise, which is expressed as 

{ci{t)cj{t + u)) — hSijS{u) . 

with 6{u) the delta function. It should be noted that the 
validity of the above form of white noise is not easy to 
justify and it is nothing else than the working hypothesis. 
The white noise is introduced to express that the random 
magnetic field is correlated on time-scales much smaller 
than the characteristic response time of the pseudo-spin 
system. It is assumed that (c^) — 2h, and its probabil- 
ity distribution may be given by the standard Gaussian 
functional form 



may be 
[ll[l4 



P[c{t)] = exp 



(23) 



Using this distribution, the propagator K between two 
definite spin states at two different times, is given by the 
functional integral: 



X[J(t)|J(0)] 



A(J(t))-c(t)] 



X exp 



2h 



dt 



V[3]V[c{t)] (24) 



with 6 being the Dirac delta- functional. We need to 
calculate the Gaussian integral with respect to cit), in 
which the functional Jacobian factor should be carried 
out. The details of this is briefly given in Appendix B 
and the result is 



K[3{t)\3m = / exp 



1 

2h 



d3 

'dt 



A(J) dt 



(25) 

Expanding the squared term inside the exponential term, 
one can put this functional integral in a familiar form of 
the path integral for a particle in the vector potential A 
and the scalar potential V = ^ ■ The parameter h just 
corresponds to the Planck constant. We can formally 
write the above functional integral by the quantum me- 
chanical path integral, that is, by using the imaginary 
time T = it 



K 



IT 

+iA V}dTWJ(t) (26) 

ar 



B. The Fokker-Planck equation 

The derivation of the Fokker-Planck is derived most 
directly by using the above path integral (HH). If we in- 
troduce the "wave function" i^{3, t), we have the integral 
equation: 

«'(J,t)= /x[J(T)|J(0)]vI/(J,0)dJ(0) (27) 



and following the standard procedure, we obtain the 
Schroedinger equation: 



OT 



-(P 

2^ 



-iA)2 + y 



with P 



a 



tTjj _ -ih\/. Returning to the real time; 
namely, i-^ = —■§['■ then writing ^ ^ P: we arrive at 
the standard form of the FP eg nation [20j: 



^ = hW'P - V • (AF) 



(28) 



This can be rewritten as the continuity equation that is 
written as = 0, where s denotes the probability 

current s = —hVP + AP and the components in polar 
coordinate (0, $) are 



se 



dP P 
h dP P 



T]J- 



sin 8 9$ 1 + 77272 
and the FP equation is written as 



dH 
dH 



1 



d 



1 dH 
sin 9 9$ 
dH 

sin e a$ 



sin Q dQ 



(29) 



We here consider several general consequences from the 
FP equation. 

(i): The stationary distribution: We consider ^ — 0, 
and we put ansatz of the Boltzmann distribution: P{3) — 
exp[— with /3 the inverse temperature /3 = -j-^- By 
substituting this into the RHS of the FP equation, we 
get 



){\/'^H ~ f3{WHf} ^ 



From this it follows the relation: 



/3h- 



1 + 77^72 







(30) 



This relation is just the fluctuation dissipation relation, 
which establishes the relation between the dissipation co- 
efficient 77 and the diffusion (fluctuation) coefficient h. 

(ii): Evolution equation for average value^^: The 
FP equation enables us to evaluate the average of func- 
tions on the sphere; F{3): that is given by (F) = 
J F{3)P{3)d3: 



d{F) 
dt 



h{V^F) - (VF • A) 



(31) 



where use is made of partial integration. Here as an in- 
teresting example, we take up F{3) = Jf{i — x,y,z) . 
Further we choose A = 7J, which represents the relax- 
ation effect, so we have for {F), 



dt 



2h-j{j!) 



(32) 
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This is solved as 



2^/, r 7 IN 

-d-exph^t]) 



(33) 



Namely, this shows a typical relaxation behavior of {Jf 
leading to the asymptotic value {Jf{oo)) — 



IV. SIMPLE EXAMPLES: THE EFFECT OF 
PINNING POTENTIAL 

The FP equation is of peculiar form that is somewhat 
different from the usual one that consists of the potential 
term, so we need to invoke specific techniques to deal with 
this. We now consider an example that can be treated 
by analytic as well as approximate ways. 

As the Hamiltonian we take the type given by the form 
such that it is given as a function of Q: That is given by 
the pinning potential ()17|) which comes from the pinning 
center located at the north pole. As a concrete form 
we consider the case that the profile p is given by the 
Gaussian [2l|: p = exp[— ^"""^^ ] with the vortex size a, 
hence we have 



H{e) = Vbexp[-a2(l -cose)/a] = VaH. 



(34) 



As another form, it is possible to adopt the potential 
arising from two pining centers; one is located at the 
north and the other is at the south pole: and is written 



H{Q) = Vo|exp[ (l-cose)]+exp[ (l + cosG)]} 

a a 

In the following we consider the two cases without or 
with the effect of dissipation. 

(i) The case that there is no dissipation: rj = 0: We 
examine this case in the semiclassical way based on the 
classical equation of motion; the "Lagrangian" in the ((25)) 
is is 

jr ^ 1 j2 [e^ + sin^ e($ + (35) 
2 J sm B dB 

We here examine a special solution such that 8 = 0, 
which leads to O = Oq. Then, the functional integral is 

K = y"exp[-^^ sin2eo($ + K)2dt]X>($) (36) 



with 



- r ^ 

K = — 
/ on 



Vo dH. 
J sin 8 dB^e^Qo 



This is equivalent to the path integral for a particle on a 
circle; the corresponding FP equation becomes 

dP d^p dP 

h^r- = — ^ ^ hK—— (37) 

dt sin2eoa$2 a$ ^ ^ 



which is a modified form of the diffusion equation. We 
assume the solution to be a form 



P{^,t) 



n— — oo 



exp[-A„t]/„($) 



(38) 



Here /„($) satisfies the eigenvalue equation 



sin^ So d<^'^ d$ 

We put /„($) — exp[m$] (n — integers), and the eigen- 
value is obtained as 



A,] 



sin^ 9o 



n + inn 



By substituting this into the (|38p . we have 



P($) 



OO 

E 

n— — OO 



exp[-( 



sin^ 



-n^ + inK)t] exp[OT$] (39) 



which turns out to be the theta functionf23|. The ex- 
pectation values for (J) = J^^ JP{^,t)d^ is calculated 
as 

ht 



(Jx) = Jsin0oexp[- 



sinGo 



cos Kt 



(Jy) = =Jsin0oexp[- 
(Jz) = Jcos9n 



ht 



sin 00 



sin Kt 



(40) 



The damping behavior of {Jx,Jy) shows is similar to the 
diffusive behavior of spin in magnetic systems plj . 

(ii) The case that there is dissipation: rj ^ 0. 

For this case, we look for the distribution function as 
a function of Q only P = Pi&) , so the stationary FP 
equation is given by 



Id, dH , 
sine-P) 



_ji <L di^ 

sin e de ^''^'^ " de 1 + 77^2 sin ede 

which is rewritten as 

1 d I . ^dP e d dH 

(smB — ) (smB P) = — AP 

sinBdB^ dB'^ sinBdB^ dB ' 

(41) 

Here we put e = i^i^^^iji-j and X = ^. This can be 
regarded as an expansion parameter. Noting that e = 
h{^+li^j'^) turns out to be e = Vq/? as a consequence of 



-EP 



the fluctuation-dissipation relation (j30p . we can treat the 
expansion scheme in terms of the expansion parameter e 
(a): We first consider the perturbation procedure: this 
may be relevant for the case that Vb/3 << 1, which means 
the temperature is higher enough the interaction strength 
satisfies Vq << k^T. Hence expanding in terms of power 
series with respect to e: 

P(B) = p(°)(B) + p(iHB) + --- 
A = Ao + Ai H 



6 



Substituting this into the above eigenvalue equation, we 
can determine each order by iteration procedure: the 
zero-th order 

jL^^„„e^).-A.P<»> (42) 
sm H dB aB 

and the first order: 



h d ^dpW 



6 d , „ dH , _/n\ 

sinBdB' dB ' sinBdB^ dO' 

= -AiP(") - AoP^i^ (43) 



The zeroth order gives the Legendre polynomiah P'^ = 
P/(B) and the corresponding eigenvalue Aq = l{l + 1) 
with I = 1 ■ ■ ■ . Let us choose a particular namely, 
= Pi(B) ; then the first order term is calculated to 

X^=eJ F,(B)-^^(sinB^)Pz(B)df7/ J P^ie)dn 

(44) 

P*^^) can be expanded in terms of Pfc(B): P^^-* = 
J2k C'fePfc(6), then we obtain the coefhcients Ck 



C, = 



Vki 



Ao(0-Ao(fc) 



(45) 



with the matrix element: 



Vki = e f p,(B)-^-^(sinB^)Pz(B)dr! 
J sm B dB 



dB- 



for I k. Using this, he second order shift of the eigen- 
value is calculated 



^kl 



k^l 



Ao(0-Ao(fc) 



(46) 



(b): We next examine the case that e ~ 1, namely, 
Vq — ksT. In this case one has to use non-perturbative 
procedure [ll|. In order to carry out this, we rewrite the 
FP equation in terms of the variable ^ = cosB: 

and if noting the relation 

r fvi d r ~ , d dH 
expHi7]-exp[ei/] = -+e— 



we have 
d 



d 



{ (1 - exp[-ei?] ^ ( exp[eH]P) } = - AP 
This can be further reduced by putting Q = exp[ei/]P, 



This equation can be translated to the variational prob- 
lem; namely, we introduce the "action" ; 



i[Q] 



(l-e')exph6i/](^)'d^ (48) 

L "S 



and the constraint coming from the normalization; 



exp[-eiJ]g^dC = N[Q] = 1 



(49) 



Taking the "Euler-Lagrange equation" 5{I[Q]-\N[Q]) = 
0, we recover the above eigenvalue equation with A is 
nothing but the Lagrangian multiplier. We see that there 
is a trivial solution Qo = constant with the zero eigen- 
value Aq = 0. Starting from this as the lowest solution, 
we have an algorithm to get the subsequent solutions; 
Qn{n = 0, • • • ) by noting the orthogonality condition; 



exp[-eiJ]gg„dC = 



(50) 



As a concrete case, we exemplify how to construct the 
second lowest state by taking a trial function given by the 
quadratic function Qi(f) = A + P^ + C^^. Substituting 
this into ([SO]) , we have the linear relation between A, B, C 
and C can be expressed by A, B. Next (|49)) gives the 
quadratic equation for A^B. Finally (|48| gives another 
quadratic form for A^B. Then, the problem is reduced 
to the minimization of the (j48p written in terms of the 
quadratic form under the constraint (1491) written by the 
quadratic equation. Having solved these, we arrive at the 
solution for Qi. This procedure may proceed to step by 
step and we can have a series of the stationary solution 
of the FP equations. 



V. EXTENSION TO NON-SPHERICAL 
VORTEX 

Here we discuss that the non-spherical vortices can be 
readily obtained as variants of the vortex on the sphere 
that has been given above; we are concerned with the 
vortex on 2-dimensional plane and the vortex on pseudo- 
sphere. 



A. The case of planer vortex 

This can be regarded as a limiting case of the spherical 
vortex. For the vortex center coordinate R = {X, Y), the 
equation of motion for R is obtained as a special case of 
the infinite radius; a oo; 



dR, dH 
mpoMkx— )^ — 



(51) 



^{(l-e')exp[-eff]^} =-Aexp[-eiJ]Q (47) 



where k is the unit vector perpendicular to (x, y) plane. 
The Langevin equation with dissipation together with 
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the white noise c is given by 
dR 



B 



(52) 



Following the same procedure as the case of spherical 
vortex, we obtain the FP equation: 



dP 



V^P - V • (BP) 



(53) 



where V 



/ a d 
plane. In the above, B is 



is just the nabla in two dimensional 



1 



, dH dH, 
,("Wkx— ) 



This is essentially same as that has been used in the 
case of superconductivity vortex • More detailed and 
complicated treatment was given for the actual super- 
conductors in connection with high superconductivity 
[25I [26| . Our treatment here provide a general way to 
clarify the stochastic behavior of planer vortex. The de- 
tails of the argument will be given elsewhere. 



B. The case of vortex on pseudo-sphere 

Next we consider a point vortex on the pseudo-sphere 
(which is denoted by PS^) which has the surface with 
constant negative curvature. The superfluid coated on 
such pseudo-sphere may be realized in actual conditions 
as a local part of the superfluid that is put on the surface 
of complicated shape and this part can be regarded as a 
pseudo-sphere. 

The point of the surface is described by the equation 
-|- 2/^ — = — a^, which is parametrically written as 

X = sinh 6 cos (t>,x ~ sinh 9 sin (j)^ z — cosh 

Let R = R/a — (sinh 9 cos $, sinh 8 sin cosh 9) be 
the vortex center, which satisfies [27| 



rfR - dH 

LLll = — R * — 

dt dH 



(54) 



Here the scalar and vector product , denoted as *, are 
defined by 

a • b = aihi + 0262 — 03^3 

a*b = (0263 - a3fe2,a3fei - a2^3,-(ai&2 - a2&i)) 
and the triple product formula 

a * (b * c) = (a • b)c — (a • c)b 



PS^ is trivial in a topological sense, namely, this is iso- 
morphic to 2-dimensional flat Euclidean space. We in- 
troduce the pseudo spin; K — Kli with K = ^V,. The 
norm of K is given by = —K^. K is not quantized, 
which is a consequence of the non-existence of topologi- 
cal invariant for PS^ . According to the same manner as 
the sphere vortex, the Langevin equation with the dissi- 
pation is obtained as 



dK 
i(K) 



A(K) 
1 



1 



dH ^dH, , , 



Here note that c, which is the random force, has the norm 
such that = cf + — Cg , which is a consequence of the 
non-Euclidean metric on the pseudo-sphere. By taking 
account of this fact and following the same procedure as 
the spherical vortex case, we arrive at the FP equation; 



dP 



sinh 9 
where the current is 



(56) 



se 



S4, 



, dP 
h dP 



P 



P 



sinh 9 9$ 1-7^2X2 



^^99- 
dH 



1 dH \ 
sinh 9 a$ j 
TjK^ dH 
sinh 9 9$ 



We can develop the approximate procedure in the same 
way as the spherical vortex. This may be realized by 
finding the pinning potential, which can be obtained by 
replacing cos 9 — > cosh 9 (fT7|) . For developing the per- 
turbation, the zeroth order (unperturbed) distribution 
fimction satisfies 9, 



1 d . dP 



AP 



and the eigenvalue is given by A = p{p + 1) with p is the 
real number and the eigenfunction is written in terms of 
the integral representation (28j : 



Pp (cosh 9) = 



1 

2^ 



(cosh 9 + sinh 9 cos S,yd^ 



Further, for the non-perturbation case, we can construct 
the similar equation with (|47| as well as the (|48l) . which 
is simply obtained by replacing ^ — cos 9 by ^ = cosh 9. 
Then we can apply the variation principle for I[Q] leading 
to the approximation solutions for the FP equation. 



VI. SUMMARY 



Note that the 3-rd component of the scalar product has 
a minus sign compared with the spherical case. Further- 
more noted is that for the case of vortex on PS'2 , there is 
no counterpart of the topological quantization, because 



We have presented a general formalism of the stochas- 
tic theory of random behavior of the quantum vortex that 
is coated on the two-dimensional sphere. The random- 
ness are caused from several origins. We have mind of 
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the conventional assumption; temperature fluctuations. 
For the case of mesoscopic size we have in mind, the 
fluctuation of temperature is more appreciable than the 
macroscopic system from the relative sense. As the other 
origin, we can take into account of the one that is intrinsi- 
cally built in the spherical geometry, for which we expect 
to have the irregular deviation from the ideal sphere in 
actual situation 29, 30J. 

To formulate the random behavior, we follow the anal- 
ogy with the spherical vortex and spin. This analogy 
seems natural; since the vortex carries angular momen- 
tum of the fluid in a form of "lump" which is the small 
region that angular momentum is concentrated. From 
the technical point of view of stochastic theory, the 
present approach is closely connected with the magnetic 
system pTij. Having noted the analogy with spin, we can 
write the Langevin equation for the Brownian motion 
of vortex, which can be converted to the Fokker-Planck 
equation using the functional integral method [H, 
Though this present approach shares its idea with the 
FP equation for magnetic system[ll|, the present func- 
tional integral provides with the more direct means to 
access to the FP equation, that is quite different from 
the stochastic integral method (Ito integral) used in [ll| . 
As a particular case, the approximate scheme has been 
developed for the case that the vortex Hamiltonian is 
given by the pinning potential. This case may be sig- 
nificant for practical analyses, since the pinning plays a 
crucial role for the vortex motion. 

Finally, we gave an extension to the case of the variants 
of the spherical vortex; the planer vortex and the vortex 
on the pseudo-sphere. The vortex on the pseudo-sphere 
may be regarded as a step towards the quantum vortex 
that is defined on the more complicated two-dimensional 
manifold. 
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Appendix A: Calculation of V x j 

We rewrite the surface integral of the right hand side 
of ([5]) by using the Stokes theorem: 



/ (V X j)(icr — / mpvds 
Js Jc 



(Al) 



with C being boundary of S . By noting the fact that if 
C is chosen such that it is sufficiently remote from the 
vortex center, p — > poj we have 



mpvds = mpo / v ■ ds — mpo / ujda 
c Jc Js 



with the vorticity a; = V x v. 

Next we give an explicit form for the velocity field. 



which is given by 

v = ^(fxV/), (A2) 
where the vortex charge p is introduced and / is 

/ = log |r - R(t)| loga + log |f - R(t)|. (A3) 

Here f = r/a,R = R/a, that is , the vector on unit 
sphere: 

R = (sin cos $, sin Q sin cos O) 

At this point we adopt the polar coordinate system: 
f = (sin 9 cos <j), sin 6 sin (p, cos 9) and the velocity field be- 
comes 



sm & oca aO 



(A4) 



Hence the vorticity is given by w = V x v = (V^/)f with 
the Laplacian on 2-sphere 



1 d , . Bf 1 



sin 6* 961^ de' sm^edch^ 



where we note log |f — R(i)| = 6{r — R(i)). From the 
above expression the vorticity for the spherical vortex 
directs to the normal of sphere. By substituting this 
into the Lagrangian equation of motion, we arrive at the 
Kirchhoff equation. 

Vortex equation on a pseudo-sphere 

Here we give a sketch for the vortex on pseudo-sphere 
PS^ 

The procedure of calculation is completely same as the 
case of S^; we need only the following changes: 

•' ^ ' de' smhO dcj)' de smhedcj)^ 

where 9, (j) means the unit vector in the direction of {6, 4>). 
Noting that the velocity field is given by replacing sine by 
sinh 9 and hence the vortex becomes w = V x v = (V^/)f 
with the Laplacian on pseudo-sphere 



1 d 



V^/=-r7I7j((sinhe^) 



1 d^f 



sinh 61 a de' sinh^ edcj)^ 

Hence, we obtain the equation of motion in the same form 
as the spherical case by noting the difference of vector 
product. 



Appendix B: Reduction of functional integral 

Here we show how to reduce the Gaussian functional 
Ltegral. We introd 
functional identity: 



integral. We introduce G = ^ -|- A(J). First we use the 



dct [^] J n ^[G(J) - c]D[3{t)] = 1 (Bl) 



9 



then we have the functional integral including the the 
Jacobian factor I3lll32ll: 



l[S[Git)-cit)]D[Jit)] 



D 



X exp[ 



c^{t)dt]D[c{t)] 



(B2) 



Here we can remove the constraint that attaches to 
the det, namely, G = c. Because, it can be obtained 
by the well known the relation for the delta function 
c{x)6{x — a) = c{a)S{x — a) and the determinant fac- 
tor is generally regarded as a functional of 3{t), so it 
does not depends on c(t) any more. Further, using an 



integral identity S[g{x)] = J exp[iXg(x)]dX^ the above is 
written as a functional integration over three variables 
J, c. A: Thus we have 

K ^ j exp[i j X{t) ■ {G{t) - c{t))dt\ 
D f 

X exp[-y / c^{t)dt]D[\{t)] 



X Aei[^]D[c{t)]D[m] 



(B3) 



By accomplishing the two Gaussian functional integrals 
for c and \{t) successively, we obtain ([25)) . 
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